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Предложена модель процесса электронно-лучевой наплавки покрытий с модифицирующими частицами с учетом физико-хими- 
ческих превращений. Процесс растворения описывается на основе формально-кинетического подхода; процессы плавления и 
кристаллизации ~ в рамках теории двухфазной зоны. Результатом численного решения задачи является фазовый и химический 
состав образующегося покрытия, размер зоны термического влияния и ванны расплава в зависимости от параметров, характе- 
ризующих режим наплавки. 


В настоящее время широкое распространение 
получили электронно-лучевые технологии обра- 
ботки материалов, нанесения покрытий, получе- 
ния неразъемных соединений и т. п. В большинстве 
случаев такие технологии являются трудоемкими и 
дорогостоящими, поэтому для предварительного 
выбора технологического режима большую роль 
приобретают теоретические исследования техно- 
логических процессов. Известные к настоящему 
времени математические модели высокотемпера- 
турных технологических процессов, в основном, 
обращают внимание на чисто тепловые процессы: 
теплопроводность, плавление, кристаллизацию, 
испарение. По значению температуры, достигав- 
шейся в зоне обработки, авторы различных работ 
судят о фазовой структуре материалов. 

В данной работе предпринята попытка учесть 
влияние реальных физико-химических процессов 
(таких как растворение твердых частиц в расплаве, 
химические превращения в твердой и жидкой фа- 
зах и др.) на фазовую структуру формирующихся 
покрытий при различных режимах термической 
обработки. 


Математическая постановка задачи о формиро- 
вании фазовой и химической структуры покрытия 
в процессе электронно-лучевой наплавки с моди- 
фицирующими частицами [1] включает уравнение 
теплопроводности [2] 
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и уравнение кинетики растворения для частиц 

^ = Я,„-Ф (Тр П р ,-У ( 2 ) 

где с, р, X - теплоемкость, плотность, коэффициент 
теплопроводности (зависящие от температуры); 
слагаемое ест( Г 4 - Г 0 4 ) описывает теплоотдачу с по- 
верхности ванны расплава по закону Больцмана; 
Т - температура; 0 А - тепловой эффект химиче- 
ской реакции; () р - тепловой эффект растворения 
частиц в расплаве; р р - объемная доля нерастворив- 
шихся частиц; Я - толщина пластины, по которой 
движется источник или эффективная ширина обла- 


сти, соответствующей зоне термического влияния 
вглубь массивного образца (во втором случае в ( 1), в 
принципе, можно учесть потери тепла вглубь образ- 
ца за счет введения эффективного коэффициента 
теплоотдачи без введения третьего измерения). 

Энергия в эффективном источнике д е , движу- 
щемся по обрабатываемой поверхности со скоро- 
стью V, распределена по закону 

_|о, у>\\ 

1с ІЯо ех Р (~(х-Ѵі) 2 /а 2 ,), у<И о, (3) 

что соответствует пилообразному режиму сканиро- 
вания [3] электронного луча и включает ширину 
сканирования й 0 и величину а„ приближенно рав- 
ную радиусу электронного луча. 

Плотность потока частиц, поступающих в рас- 
плав, удовлетворяет уравнению 

Ч,„ = «Яо ехр (-((х-х а -Ѵ() 2 + у 2 ) /а 2 ), (4) 

где д и0 - максимальная плотность потока частиц; 
величина а р определяется радиусом трубки, через 
которую подаются частицы, х а - расстояние от 
максимума энерговыделения. 

Кинетическая функция ф, характеризующая 
скорость растворения, в общем случае зависит от 
температуры, дисперсности, растворимости ком- 
понентов друг в друге и локальной гидродинамиче- 
ской обстановки [4, 5]. В модели принята зависи- 
мость ф(Т)=ф р (гі р )к^ехр(-Е а /КТ), где Я - универ- 
сальная газовая постоянная, к 0 - предэкспонен- 
циальный множитель, Е а - энергия активации. Со- 
гласно литературным данным [4, 5], если скорость 
растворения твердых частиц в расплаве лимитиру- 
ется процессом диффузии, то Е а близка к энергии 
активации диффузии. Если перенос элементов из 
твердой фазы в жидкую сопровождается химиче- 
ским взаимодействием элементов, то Е а является 
некоторой эффективной величиной, зависящей от 
преобладания тех или иных процессов. Вид кине- 
тической функции зависит от процессов, опреде- 
ляющих скорость растворения на микроуровне. 
Константа к„, в соответствии с известными моделя- 
ми растворения, определяется характеристиками 
локального гидродинамического течения в окрест- 
ности частицы, попавшей в расплав. 
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Технические науки 


Если частицы нерастворимы в расплаве, то <д=0 
(или С с =0); если вещества, из которых состоят ча- 
стицы и основа, бесконечно растворимы друг в 
друге или химический состав образующегося ра- 
створа не представляет интереса, то (4=0. Вариан- 
ты модели электронно-лучевой обработки, соот- 
ветствующие нерастворимым частицам и неогра- 
ниченной растворимости материалов частиц и ос- 
новы друг в друге, проанализированы в [6-9] . Если 
нельзя пренебречь химическими реакциями в ра- 
створе, то уравнения (1-4) должны быть дополне- 
ны кинетическими уравнениями, а уравнение те- 
плопроводности - объемным тепловыделением 
вследствие химических реакций 

Ол =рІА(<ЧМ), (5) 

к = 1 

где Іі к - парциальная энтальпия (образования) к-то 
вещества, р - плотность материала основы, г/ к - 
массовые концентрации индивидуальных веществ 
и соединений в смеси, п - количество веществ, 
участвующих в реакциях. 

Теплоемкость в (1) вычисляется по различным 
формулам, зависящим, как и параметры растворе- 
ния, от типа изучаемой системы [2]. 

В качестве модельной системы с растворимыми 
частицами и химическим взаимодействием выбе- 
рем систему А1+Си. Процесс растворения этих ме- 
таллов друг в друге сопровождается химическими и 
фазовыми превращениями. 

В области концентраций 0.. .0,7 вес. % меди на 
диаграмме состояния [8]) мы можем выделить 2 хи- 
мических соединения: СиАІ и СиА1„ образование 
которых идет в соответствии с реакциями 
Си + А1 -» СиАІ, СиАІ + А1 -» СиАІ,, 

Си +2А1 — » СиАІ,. (6) 

Следовательно, запишем кинетическое уравне- 
ние для каждой реакции 

5ц, 5ц, 5ц 3 ... 

= -ф 1 -(Рз,-^ = ф 1 -ф,, — = ф 2 +ф 3 , (7) 

где фі - функции, зависящие от температуры и кон- 
центрации веществ, участвующих в реакциях; сг - 
скорость изменения концентрации меди в растворе: 

_ 5С,„ _ Ф(Ц„Г)М А ,М Си 
5/ (^М Си + (1-^)М А1 ) 2 ’ 

где М Си , М м - молярные массы меди и алюминия, 
2; - объемная доля материала частиц в растворе, 
удовлетворяющая уравнению 

§ = ИГ 


С„, - массовая доля меди в растворе, следующая из 
соотношения: 


С = 




<■> м Са + (і - )М А 


( 8 ) 


В соответствии с законом действующих масс и 
законом Аррениуса, запишем 
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Константы скоростей реакций к ш , значения 
энергий активации Е ш и те плот реакций 0, (/=1,2,3) 
определим с помощью известных термодинамиче- 
ских формул. Так, константу скорости реакции 
можно найти с помощью формулы [9] 
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где Л^=6,022Т0 23 - число Авогадро, молы 1 ; 
й=6, 626176-КТ 34 - постоянная Планка, Дж-с; 
А5" еак - энтропия реакции, Дж/(моль-К). 

Для каждой из рассматриваемых реакций имеем 

4 = ^АІ + *8*Си — ^АІСи > *^2 = ^АІСи + *^Си — ^СиАІ 2 > 

*^3 = *^2А1 + *8си ~ *^СиА1 2 ■ 

Из справочных данных [10] удается определить 
все значения энтропии веществ и соединений, кро- 
ме энтропии фазы СиАІ,: 5^=28,33, 5 Си =33,24, 
5 С иАі = 32,72 Дж/(моль-К). Отсутствующее значение 
энтропии образующегося соединения рассчитыва- 
ем по формуле Истмена [11, 12] 

8 ж = К (| 1п 4 Р + 1п К Р - § 1п 4л ) + а, 

где Д. р - средний атомный вес, (т. е. молекулярный 
вес, деленный на число атомов в молекуле) 
(Адіаі 2 — 65,02- 10~ 24 ); Ѵ ср _ средний атомный объем 
( ^сиАі 2 = И ,0872Т0~ 30 ); Т ш - температура плавления сое- 
динения; а - константа, равная 12,5+2 кал/(град-г-атом). 
Расчеты дают 5 , СцА] =-1,223 Дж/(моль-К). 

В результате имеем 

А: 01 = 1,64-10 15 ; А: 02 = 6,8-10 16 ; к ю = 2, МО 18 1/с. 

Энергию активации каждой реакции можно 
найти из формулы [9] 

Е ш =Т(8 і+ К). 

Следовательно, имеем: 74=33406, 74=61249, 
Е аі = 8 7 1 8 5 Дж/ моль. 

Теплоту образования (энтальпию) для чистых 
веществ А1 и Си принимаем равной нулю, т. е. й Си =0 
и А а 1 =0, а энтальпии соединений - к см = 39,81, 
Асыи. — 39,75 Дж/моль. 

Подставляя (7) в (5) и группируя слагаемые с 
одинаковыми (р,, получим выражения для теплот 
реакций Дж/моль (/=1,2,3): 

а=^(Лс иА1 -Л Сц )=5,6Ы0 3 , 

1Ѵ1 Си 

02={^(Ьсш 2 -Ьсш) = - 11, 22 - Ю 3 , 

1Ѵ1 Си 

а=^(4иА, 2 -4и) = -5,6М0 3 . 
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Температуры ликвидуса Т ІЩ и солидуса ^*>1 ап- 
проксимируем линейными зависимостями на нес- 
кольких интервалах [8] 


Т . = а +Ъ -С ; 

80 І 8 X т ’ 

Т\щ = а , +Ц • с„ . 

(9) 

Граничные условия имеют вид: 


х — > ±оо : 

^ = о- 

дх ’ 

(Ю) 

0: 

II 

О 

(П) 

у^со: 

^=0 

ду 

(12) 

и отражают отсутствие источников тепла при 

бес- 


конечном удалении от ванны расплава и на торце 
пластины, принятом за начало отсчета. 

В начальный момент времени имеем условия 
( = 0: т = Т о ,г] 1 =0, г\ р = 0 ,?^ = 1 . ( 13 ) 

Задачу (1-4), (7), (9-13) решали численно по 
неявной разностной схеме второго порядка ап- 
проксимации по пространственным шагам и пер- 
вого по времени с использованием расщепления по 
координатам и линейной прогонки. 

Анализ результатов численного исследования 
задачи показывает, что по мере продвижения ис- 
точника вдоль обрабатываемой поверхности уста- 
навливается квазистационарный режим. Этот ре- 
жим имеет место и для нерастворимых частиц \У 
[13], и для растворимых частиц № [2, 6, 7], и для 
случая ограниченной растворимости частиц Си, 
рис. 1. 

Максимальная температура Т тш вдоль оси дви- 
жения источника, размер зоны термического влия- 
ния у Т , х т и ванны расплава у т , х т ; доля частиц на 
границе зоны термического влияния за движущим- 
ся источником 77 рй (в точке А, рис. 1, б) по истече- 
нии некоторого времени практически перестают 
изменяться. Реально расчет идет до тех пор, пока 
квазистационарный режим не установится с задан- 
ной степенью точности. Как в случае растворимых 
частиц, так и в случае нерастворимых частиц, вели- 
чины у Т , х т , у ю х я , Г тах , т] ря зависят от технологиче- 
ских параметров. 

Для оценки свойств образующегося покрытия 
наибольший интерес представляет зависимость до- 
ли частиц в покрытии от параметров источника для 
заданных свойств материалов. Распределение п Р 
для различных Г и различных значений плотности 
мощности источника % показано на рис. 2, а. Оче- 
видно, что чем выше плотность мощности источ- 
ника, тем больше частиц растворяется в ванне рас- 
плава. Рис. 2, б, в, дают представление о распреде- 
лении химических соединений вдоль оси движения 
электронного луча при различных значениях тех- 
нологических параметров. Состав (доля нераство- 
рившихся частиц и доли соединений в растворе) на 
квазистационарной стадии принимаются за состав 
получаемого покрытия. 
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Рис. 1 . Процесс электронно-лучевой наплавки: а) распреде- 
ление температуры вдоль оси движения источника; 
6) форма и размер ванны расплава и зоны термиче- 
ского влияния. До =5- ІО 3 Вт/см 2 ; д т0 =7,5 1/с; х а =0 

На рис. 3 показано влияние плотности мощно- 
сти источника и плотности потока частиц на фор- 
му и размер ванны расплава и зоны термического 
влияния для момента времени, когда процесс мож- 
но считать установившимся. 

Ванна расплава закрашена темно серым цветом. 
Ширине зоны термического влияния, наблюдаю- 
щейся экспериментально, соответствует макси- 
мальный размер области светло-серого цвета вдоль 
оси у, определенный по температуре 700 К. 

Увеличение плотности мощности источника су- 
щественно изменяет картину. Увеличение же плотно- 
сти потока частиц практически не оказывает влияния 
на качественную картину развития процесса. 

На следующем рисунке (рис. 4) представлена 
зависимость объемной доли образующихся соеди- 
нений в точке максимальной температуры по оси х 
от параметров, характеризующих режим наплавки. 
Зависимость состава покрытия от скорости движе- 
ния внешнего источника вдоль обрабатываемой 
поверхности аналогична: чем выше V, тем ниже 
температура расплава, тем с меньшей вероятно- 
стью образуются соединения СиАІ,, СиАІ. 

В результате в процессе электронно-лучевой 
наплавки образуется либо композиционное по- 
крытые с некоторой объемной долей частиц (зави- 
сящей от параметров V, д 0 , д т0 ), внедренных в ма- 
трицу, химический состав которой зависит от тех 
же параметров, либо практически гомогенное по- 
крытие, если частицы растворяются полностью. 
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Рис. 2. Распределение: а) объемной доли частиц, б) образую- 
щегося соединения, в) материала основы в различные 
моменты времени і, с: 1) 0,5, 2) 3,2, 3) 4,8, 4) 7 для раз- 
ных значений плотности мощности источника д 0 . 
Ѵ=1 см/с; Дм =1,5 1/с;х а =0; г/ 3 - СиАІ 2 ; г\ь~ АІ; р р ~ Си 
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Рис. 3. Влияние плотности мощности источника и плотности 
потока частиц на форму и размер ванны расплава и 
зоны термического влияния. Указанные моменты 
времени соответствуют времени установления ква- 
зистационарного режима или превышают его. 
а) д 0 =3, 5-107 Вт/см 2 ; д^=1,8 1/с; б) р 0 =51О 3 ; д„, 0 =7,8; 
в) До = 7' ІО 3 ; д т /о = 1,8 



Рис. 4. Зависимость объемной доли соединений, образующихся в растворе, и доли частиц в покрытии от плотности мощности ис- 
точника До и плотности потока частиц д^. Прямые линии д 3 ~СиАІ 2 , пунктирные ~ д 4 - АІ ; штрихпунктирные ~ р р ~Си. а)д т0 : 
1-0,9; 2~ 1,5; 3~1,8 1/с; б) д 0 : 1 ~ 3,5 ■ Ю>; 2-5,0 ■ 1(Г; 3~7,0 ■ 10> Вт/см 2 


123 




Известия Томского политехнического университета. 2006. Т. 309. № 6 


Таким образом, в работе предложена модель 
технологического процесса электронно-лучевой 
наплавки покрытия с модифицирующими части- 
цами для модельной системы с образованием хи- 
мических соединений. Впервые теоретически уста- 
новлена связь фазового и химического состава по- 
крытия с технологическими параметрами, характе- 
ризующими режим наплавки. Подобные зависи- 
мости наблюдаются и в экспериментальных иссле- 
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13], может служить основной для формулировки и 
решения задачи оптимизации технологии наплав- 
ки [1] с целью получения покрытия требуемого фа- 
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заданные физико-механические свойства. 
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